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Abstract 



It is shown that the self inductance of a wire loop may be written as a curve 
integral akin to the Neumann formula for the mutual inductance of two wire loops. 
The only difference is that contributions where the two integration variables get too 
close to each other must be excluded from the curve integral and evaluated in detail. 
The contributions of these excluded segments depend on the distribution of the 
current in the cross section of the wire. They add to a simple constant proportional 
to the wire length. The error of the new expression is of first order in the wire radius 
if there are sharp corners and of second order in the wire radius for smooth wire 
loops. 

1 Introduction 

Electrical inductance plays a crucial role in power plants, transformers and electronic 
devices. The coefficients of self and mutual inductance required to quantitatively describe 
inductance belong to the field of magnetostatics. Calculating inductance coefficients with 
analytic techniques is, however, impossible except in simple cases. The mathematical 
reason for the difficulty is that the Laplace equation allows analytic solutions only for 
some symmetric constellations. There thus are only a few closed-form expressions for 
these coefficients. In practice one often is forced to use approximations, finite element 
methods or other numerical techniques. The situation is somewhat less gloomy in systems 
where the current flows in thin wires. This situation is analogous to an electrostatic 
system of point charges, where electric fleld and electrostatic energy directly follow from 
the given charge distribution, while in a generic system charge and current distributions 
also are unknown at the outset. 

The purpose of this article is to derive a new expression for the self inductance of a 
wire loop, giving self inductance as a curve integral similar to the Neumann formula for 
mutual inductance. The starting point is the expression 



for the magnetic fleld energy of a system with current density j (x), where /iq is the 
magnetic constant. [1] This expression essentially was already given by Neumann in 1845. [2] 





It resembles the expression for gravitational or electrostatic potential energy, the only new 
ingredient is the scalar product between the current elements. 

For a current density j (x) = ^ Imjm (x) corresponding to separate current loops 
with currents and normalized current densities it follows 
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If the currents flow in thin wires, then the integrals become curve integrals, and one 
immediately reads off the Neumann expression for mutual inductance of two (filamentary) 
current loops [2] 

_ Ho f dxi- dx2 
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It is plausible that there exists a similar expression for the self inductance of a wire loop, 
but we were not able to find any hint in the literature. Formally one might read off from 
equation (|2| an expression similar to equation ([3]), where the two closed curves coincide. 
But this cannot be correct, because |x — x'| now vanishes and the integral isn't defined. 
Instead we will proof 

where a denotes the wire radius and / the length of the wire. The variable s measures the 
length along the wire axis. The constant Y depends on the distribution of the current in 
the cross section of the wire: F = if the current flows in the wire surface, Y = 1/A when 
the current is homogeneous across the wire. The ellipses represents terms like O (/ioo) and 
O [hqO? /I), which are negligible for / a. 

In hindsight it is completely natural to use a cutoff of order a in the curve integral. In 
fact, the exact value of this cutoff is arbitrary, because the contribution proportional to 
lY also depends on this cutoff. The simplest way to determine Y would be to compare 
the expression with the self inductance of a long rectangle. 



2 Simple derivation 

Consider equation ^ with = 1 for a thin wire with circular cross section, radius a and 
length I. Let s denote the length along the axis of the wire. The planes perpendicular to 
the wire axis then define a projection from the bulk of the wire onto the axis, x — )■ s (x) (we 
won't consider here cases where the centre of wire curvature is within the wire, and this 
projection isn't unique any more). Then, selecting a length scale h satisfying a b I 

allows to write L = (hq/Ah) (^L + with 

L = ([ '-P^d^xd^A , (5) 
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Figure 1: A section of a wire with radius a, with a segment of length 2b, and a plane 
perpendicular to the wire axis at the center of the segment. 

The second part contains contributions from all point pairs {x, x'} with a distance along 
the axis smaller than b, the first the complement of this set (s is a cyclic quantity). For 
given X the planes at s (x) ± b delimitate the points x' contributing to the first or second 
integral, see figure L now approximately becomes a curve integral and L essentially 
consists of cylinders of length 2b. 

The strategy then is to replace L with the curve integral and to explicitly evaluate 
the contribution of the cylinders in L. The cylinders are long in comparison to the radius 
because of a <^ 6 and straight (at least most of them) because of 6 -C /. Actually the only 
requirement for the lengths is a ^ /, the length b = y/al then satisfies a ^ 6 ^ /. The 
approximation thus is exact in the limit a <^ I except in special cases. Inserting Lq for a 



straight segment from equation (A. 4) in the appendix thus leads to 



/^o^/^x.rfx'^ +f^L(^I]+Y]+... (6) 
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This expression cannot depend on the (more or less) arbitrary length scale b. The curve 
integral thus is L{b) = const — ^ln(26/6o)- But b is the only "short" length scale in 
the curve integral, and L(b) thus also is valid for 6 = a/2. The expression ^ therefore 
doesn't change if one formally sets b = a/2. Equation ^ now agrees with equation Q, 
but some questions remain. 

First of all, how accurate is formula Q? The curve integral is a purely geometric 
quantity with dimension "length" and order of magnitude I. Plausible expressions for 
the order of magnitude of the relative error are a/l, (a/l)'^ and {a/Rf, with R a typical 
curvature radius of the wire loop. Errors of this order normally are negligible (and also 
occur in the Neumann formula for mutual inductance). But the derivation of formula (|4]) 
is not as straightforward, and so what are the actual limits or exceptions? 



3 Examples and comparison with exact self induc- 
tance 

To get an impression of the accuracy we have compared self-inductances calculated with 
the curve integral with the result of a numeric evaluation of the nominally 6-dimensional 
integral in equation ([2]). This integral becomes 4-dimensional if the currents flow in the 
wire surface (skin effect, F = 0), and the results below correspond to this skin effect case. 
The order of magnitudes of the error terms identified here are corroborated below in a 
more detailed derivation of formula (141). 
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Straight segment 



The first example is a straight segment with length c and complete skin effect. This of 
course isn't a closed circuit, but it might be an edge of a rectangle. The orthogonal edges 
of the rectangle don't interact with the segment because of the scalar product j ■ j'. What 
is missing for a rectangle are the interaction terms with the opposite edges (and the small 
contributions from the corners). In this case the volume integral ^ may even be evaluated 
analytically with the result 
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while the curve integral Q leads to 
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The difference is of order O (/ioo), much smaller than /iqc for c ^ a. 



(7) 



Circular loop 

The next example is a ring with radius R. The curve integral ^ gives 

= f^oR (In {8R/a) -2 + Y)+fioO [a^/R) . 

This expression also may be found in the literature, derived with the help of elliptic 
functions and some approximations in a much more complicated way. The table displays 
the ratio of the exact inductance and L^. for some values of i?/a. 



R/a 


12 3 5 


10 




L/L, 


7.45 1.224 1.075 1.021 


1.0046 


between 1 and 1 + 0.5 {a/RY ■ 



The expression Lc is more accurate than one might expect. It gives a reasonable approxi- 
mation already for R = 3a, and the error roughly decays like O (a^/i?^). 



Rectangle 

This case is more complicated because in principle also the shape of the corners comes into 
play (curvature radius?). But the simplest thing to do is to evaluate the curve integral 
^ for a rectangle with edges of length c and d. Orthogonal edges decouple because of 
the scalar product j ■ j', and the first contribution are the terms ([T]) for the four edges by 
themselves. The second contribution are the parts of the curve integral ^ with x on one 
edge and x' on the opposite one. The condition |s — s'| > 6 is irrelevant for these cross 
terms, and one easily obtains for parallel edges of length c and distance d 

(c, d) = ^( AVc^ + d^ -Ad -Ac asinh (c/d)) , 
47r V / 

and the sum together with the y-term of equation ^ is 

L, = ^{c\n- + d\n—-{c + d){2-Y) (8) 
vr a a 



+2a/ c^ + — c asinh {c/d) — ci asinh (rf/c) + a }. 
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This expression also may be found in the hterature, with sometimes a factor 2 at the 
a-term. [3] The table displays the ratio of the numerically evaluated self-inductance L of a 
square with border length c and corners with curvature radius a and the curve integral L^. 
for different border length c, 



c/a 


5 


10 


20 


40 


80 


160 


L/L, 

(L - //io 


1.1028 

0.306 

1.0055 


1.0284 

0.285 

1.0059 


1.0096 
0.2756 
1.0017 


1.00363 

0.2706 

1.00121 


1.00145 

0.2673 

1.00052 


1.00061 

0.2673 

1.00022 



The curvature radius a is minimal in that the centre of curvature lies on the inner border 
of the wire. It is remarkable that the absolute error nearly remains constant. The last row 
of the table displays the ratio of the exact self inductance and the exact curve integral (|4]) 
(with round corners), also evaluated numerically. This expression is a better approximation 
for small c/a, where the square with round corners degenerates to a ring. 



Equilateral triangle 



The curve integral (|6j) for an equilateral triangle with edge length c consists of three times 
the expression Q for the edges by themselves and three times the interaction energy 
Lc (c, c, 120) of adjacent edges (with s on one edge and s' on the other, see appendix D 
There is no such interaction for rectangles because of the scalar product). 



= ^3c \ In ( - 



1 - In 



The table displays the ratio of the exact self-inductance L of an equilateral triangle with 
border length c and corners with curvature radius a and the curve integral Lc for different 
border length c. 



c/a 


5 


10 


20 


40 


80 


160 


L/L, 

{L - Lc) //io 


2.77 
0.8605 


1.20 
0.863 


1.057 
0.8646 


1.020 
0.8675 


1.0076 
0.8692 


1.0031 
0.884 



The absolute error (L — Lc) / is nearly constant also here. 



Parallel wires 

For a loop consisting of arbitrary long parallel wires the condition a <^ / is perfectly 
fulfilled and the curve integral gives the exact self inductance even for minimal distance 
d = 2a, 

Lc=^(\n^- + Y 



IT \ a 

This expression is the limiting case of the expression ([8j) for a long rectangle. The point 
is, that the contribution of the corners becomes negligible for a long rectangle, and that 
the replacement of the (circular symmetric) current by a line current doesn't change the 
magnetic field according to Ampere's law. Of course, the assumption of a circular sym- 
metric current distribution gets wrong in the skin effect case if the wires are close to each 
other because of additional screening currents. 

To summarize, formula ^ is rather accurate even for circuits with a linear extension 
as small as 20 times the wire radius, even if the circuit contains sharp corners. 
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Figure 2: A wire loop and distances relevant for calculating the energy of a disc at s = 
relative to other points in the wire. 



4 Error estimation 

According to equation (5) the self inductance may be written as L = (yUo/47r) (^L + L 

where L contains the short segments and L the complement. The self inductance L of 
course doesn't depend on the arbitrary segment length b. 

Let us now introduce some notation. We use a coordinate system {s, r, (f)} in the wire 
where the length s along the wire axis is cyclic with period I, and the coordinates {r, 0} 
describe planes perpendicular to the wire axis. The intersections of the planes and the 
wire are assumed to be circular, ^ r ^ a and ^ ^ 27r. The volume element reads 
dV = {1 + r cos (j)/R{s)) rdrd(j)ds, where R{s) denotes the curvature radius of the wire, 
and = at the outer border of the wire (this is possible at least locally). The volume 
element also may be written as dV = dsdA with dA = (1 + r cos (j)/R{s)) d A and area 
element dA = rdrdcf). The coordinates become cylindrical coordinates for straight wire 
segments, i.e. for R = oo. 

The current density j is normalized, that is J dA |j| = / dA |j| = 1. We will also need 
the radial moment 

a2 = (r^) = j dAr"^ |j| 

of the current distribution. In the skin effect case of course 02 = a?. 
One quantity of interest then is 

1(0) = I ds'e{\s'\-h) [ dAdA' , , ^ ^^'J'' ^""'^ ^ , , , (9) 
^ W ^' ' V |x(O,r,0)-x(s',r', </.')!' ^ ^ 

the energy of the current in the plane at s = with respect to the current at > b, see 
figure (^. The symbol 9 denotes the step function. To obtain L from L(0) requires to 
move the "reference point" s = around the wire, that is to integrate over s. The curve 
integral 

is an approximation for L (0). 

Similarly we write the short segment around the plane at s = as 

L (0) = / ds'e ib - \s'\) [ dAdA' ^^"'^t^"^^"''''^^, , (11) 
^ W ^ ' 'V |x(O,r,0)-x(s',r',0')r 
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These definitions allow to write 



L(0) = L(0) + L(0)= (12) 



s=0 



where we have added and subtracted the curve integral (0), the segment integral Lq (0) 



for a straight segment from equation (A.4) and the approximation (0) = 2 (In {2b/a) + Y) 
for Lo (0). 

The first bracket in equation (12) now is formula (|6]), the second bracket contains the 
difference of the volume and the curve integral (for |s'| > b), the third bracket is a simple 
power series in a/b, and the fourth bracket is the difference of the actual segment integral 



(11 ) and the segment integral for a straight segment. It is evident that the error becomes 



small in suitable limits, and we now want to determine the order of magnitude of the error. 
Current loops consisting of straight segments 

We first consider a current loop consisting of straight segments, each with a length compa- 
rable with the total length /, and a finite number of corners, each with a length comparable 
with the wire radius a. We show that the error is of order a. The procedure is simple in 
principle but circumstantial. 



The first error term of equation (12), L (0) — (0), contains a volume current in one 
direction and an opposite current of same strength on the wire axis. The reference point 
s = mostly will lie on a straight segment, a distance c > b away from the end of the 
segment. In the notation of appendix |A] the contribution of the segment becomes 



AL (0) = j\s'dAdA' (^7^^ - \)j j (^) j {r') = J dAdA' A, (^^^ 
1 1 \ {N^) [I 1 \ 3 (A^^) 



c2 52/4 Vc4 64/ 32 



The functions Ai and also are defined in appendix |a| here we only need (^N^^ 



An inspection of equation (A.4) shows that the second error term Lq — L^ of equation (12) 
exactly cancels the b terms in AL (0) above (the additional factor 2 in Lq — L-^ is needed 
for the same purpose on the other side of the reference point). This cancellation isn't by 
accident, it comes about because L and Lq from first and second error term are adjacent 
integrals of the same type. 

To get the actual error requires to integrate over the reference point, 

J^, V 32c4 ) Ab 32b^ ^ ^ \ / J v ; 

Here we now have selected a length b of order O (a) and neglected a'^/l against a. 

The other straight segments of the current loop also contribute to the first error term 
L (0) — (0), but are away at least a distance c > b from the reference point. It is 
shown in appendix [B] that straight segments with opposite currents in the volume and in 
the axis of the wire are equivalent to multipoles at both ends with a potential of order 
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c? I |x|^. Integrating over the reference point thus completely analogously leads to an error 
AL = O [a). It is even better, the corners (also with opposite currents) also are equivalent 
to multipoles of the same type and order of magnitude (this can be verified for instance 
from the dipole term in the next section). 

One thus may combine the ends of the straight segments and the corners to one multi- 



pole and claim that the error (a) in equation (\12n is caused by the corners. 

It is essential for this argument, that the integrals over the reference point converge, and 
the corners thus add independent error contributions. 

What if s = comes to lie on a corner? This only occurs for a lenght of order a in the 
s-integral and only contributes another error of the same order of magnitude. 



It remains to examine the last error term, the fourth bracket of equation (12), involving 
segments of length 2b = O (a). This term differs from only if s is within a corner. The 
s integral thus yields something like O (a) Lq = O (a). 

To summarize, the leading error of formula ^ in the straight segment case is of order 
/ioa, much smaller than the self inductance of order /io/ln (//a) even for moderately large 
I /a. The error grows with the number of corners and depends on the shape of the corners, 
but doesn't depend on the size of the loop. This perfectly agrees with the numerical results 
for the square and the triangle above. The leading error presumably simply is the sum of 
contributions from the corners, each a functional of order O (fioa) of its shape. In practice 
it normally suffices to know the order of magnitude of the error. 

Smooth current loops 

Another special case are smooth current loops with a minimal curvature radius Rq of the 
same order of magnitude as the loop length /. The simplest example is the circular loop. 
There are no corners in this case, which were identified in the previous section as error 
source for loops consisting of straight segments and corners. So one may expect a smaller 
error here. To identify the order of magnitude of the error for smooth loops requires more 
effort. 

The coordinates x and x' in the integral ^ above may be expanded like x = x (s, 0, 0) + 

Xi (s, r, 0), and thus 

x-x' = x^^^/ + xi - x'l, 
(x - x')^ = xl ,., + 2nXs,s' + 

with Xs g/ = X (s, 0, 0) — x(s',0,0) the distance of the projections onto the axis and |xi| 
and \x.[\ of order O (a). The abbreviations are 

= (xl-x;)^ 

The procedure now is to use the multipole expansion 
1 1 /i 1 



-'y')-7^ (S/i^^ - 3/iz/2) (13) 



+ (35/x^ - 30/1^2 + 3i/^) + O (a'^x-^,) . 
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This expansion converges for Xs,s' = l^i— xi|, the coefficients are the coefficients of the 



Legendre polynomials. Inserting the leading monopole term l/xs,s' of equation (13) into 
equation (oj) reproduces the curve integral L^(0). The higher multipole terms describe 
the difference between the volume integral L (0) and the curve integral. The s'-integral 
in the multipole terms converges and the difference thus again is mainly a local quantity. 
The nominal order of magnitude of the multipole terms L„ is a^/b^, n ^ 1. With a 
length b = y/oR the order of magnitude becomes (a/i?)"''^, and the expansion up to the 
hexadecupole {n = 4) is needed to get the {a/Rf approximation. 

It is expedient to use a coordinate-free notation for the factor (1 + rcos(f)/R (s)) in the 
volume element, 

(1 + r cos (p/R (s)) rdrdcpds = (l + Xi ■ R (s) / R^ (s)) rdrdcpds, 

where R (s) is the local curvature radius vector. Integrals over and 0' then may be 
evaluated with the help of 

(xi) = 0, 

where Pm,n is the projection operator projecting onto the plane perpendicular to the wire 
axis. This implies P ■ R = R. 

Inserting now the multipole expansion (13) into the volume integral ^ generates an 



expansion of the difference of the volume and the curve integral in the region > b. 

Dipole A contribution to L (0) only comes from the combination of the Xi from /i with 
the xi from the volume element, and vice versa for x'^^. 



cos a 



The angle a denotes the angle between the current at s and the current at s'. The integral 
over s' converges for large s' because of xo,s' ~ s' and instead of the nominal a/b there 
at most remains a^/ [bR). An interesting special case is a segment with R = oo most of 
the time and R = O (a) along a length s' also of order 0{a). The s' integration then 
reproduces the contribution Li (0) = O (a^/ |x|^) from a multipole at the position of a 
corner used above. 

For a smooth curve near s' = approximately xq.s' ■ R = Xs'_o ■ R' = s'/2. The integral 
therefore is logarithmic and not local any more. This means that the error depends on the 
global shape of the curve. But the logarithmic factor changes only slowly and the error 
gets the order of magnitude {a? /R^) In {R/b) ~ {a^ / R^) In {R/a). 

Quadrupole After integration over the angles and the radial coordinates r there 
remains 

3 / ^ ^2 /^i ^ t\ cos a 



L2 (0) = as ds'e {\s'\ - b) (2 - (n • Xo,y)' - (n' ■ - 



where n denotes a unit vector in the direction of the wire axis and the expression P = 1— nn 
for the projection operator was used. In L2 (0) one may recognize 1 — (n • xo,s')" = siia^ ip 
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and 1 — (n' ■ xo,s')^ = sin^?/^', where ip denotes the angle between the distance vector xo,s/ 
and the wire axis. With ip' = ip = 2a = s'/2R for a smooth current loop there remains 
an error {o? / R^)hi{R/h). The — 02/^^ term (from the —1) is peculiar. With the choice 
h = a/ aR it would be of order a/R, but it gets cancelled by the third error bracket Lq — 



of equation (12), see also equation (A. 2) and (A. 4) 



Oktupole The oktupole contributes at most a term of order a^/b'^. But the cos0 and 
sinip from the odd power of fi make the actual contribution smaller than the expected 



Hexadecupole The hexadecupole term is of order a'^/b^ ~ / R^ and its leading part 
must be kept. The /i factors contain a factor s\.nip ^ b/R and may be dropped. The u'^ 
leads to 

I4 (0) = ^ / ds'e {\s'\- b) [2a4 + a^3 + cos' a)] 



0,s' 



The factors cos a may be replaced with 1 because of a = s' / R. The integral converges and 
contributes an error like a'/ R^. But this ter m als o gets cancelled by the third error bracket 
Lq — Ly of equation (12), see also equation (A. 4). To be more exphcit, the Ai (N/b) gives 
an N"^ which gives 204 



4a2 after averaging over the angle 



For smooth current loops it finally follows 

An 



L= (^L^ + L^^ + O {{a^/R) In {R/b)) + 



(14) 



where the first bracket on the r.h.s is formula (|6j), evaluated with b = y aR. The remaining 
segment integral L for a segment with curvature radius R is evaluated in the appendix, the 



result (C.l) contributes another logarithmic error {a? / R^)\n{R/b), and an error b"^ / R^ 



But the latter term only cancels the 6-dependence of the curve integral L^, evaluated as 



(C.2) also in appendix C As was argued in the simple derivation of the formula above the 



cancellation comes about, because b is the lower limit in the curve integral and the upper 
limit in the volume integral. 

To recapitulate, for segment length b = y/aR there remain formula ^ and small terms 
of order O {{a'^ / R) In (i?/a))Q We now drop these small terms, since they only contribute 
to the expected error. There is no small lenght scale in th e for mula L-y + L-y and its 
6-dependence is under control for all b < \/ a R an d given by ( |C.2 ). It contains a "large" 
\? jR? ~ ajR term, which gets cancelled by (C.l). In total there only remains + 
without its "large" a/i?-term. But this is identical with + evaluated for b = a/2, 



apart from another contribution to the a^ /R^ error, see (C.2). In other words, the simple 
formula ^ approximates self inductance with an error of order (a'^/R) In (R/b) or a'^/R. 



5 Conclusions 

The curve integral ^ for the self inductance of a wire loop is only a little bit more 
complicated than the Neumann formula for the mutual inductance of two wire loops. 



^The error terms may be evaluated analytically to order /E? for a cicular loop. The result agrees 
with the result of Rayleigh and Niven in the constant current case (see ref. ^ ) and perfectly agrees with 
the numeric results above in the skin effect case. 
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The exact expression for self inductance is a 6-dimensional integral with a logarithmic 
divergence and several length scales. Nevertheless clear statements follow for the accuracy 
of formula Q, for loops consisting of straight segments as well as for smooth loops. The 
error originates from the curved parts of the loop, and is of order fiQa or ^qo^/I, negligible 
for most practical purposes. The leading error presumably even might be given as a sum 
over the corners for loops consisting of straight segments or as additional curve integrals for 
smooth current loops (see the dipole and quadrupole contributions above). The techniques 
used for error estimation also might be used for the Neumann formula. 

The self inductance curve integral can be evaluated analytically in many cases, for 
instance for current loops consisting of coplanar straight segments. But also the numer- 
ical evaluation of a two-dimensional integral with a computer program is a breeze with 
appropriate numerical libraries. 

Current distributions which are not circular symmetric also lead to formula (|4]), with 
a cutoff and a constant Y depending on the current distribution. An example are circuits 
consisting of coplaner flat strips of width w. The self inductance of such circuits is given 
in the accuracy described above by formula Mb with a = w and Y = 3/2. 



A Contribution from straight segments of length 26 



The contribution L to the self inductance in equation ([5j) is due to the interaction of 
the current in the plane s with the current in all planes s' with |s' — s| < b. This value 
depends on the current distribution in the wire and on the wire geometry within the 
segment [s — b, s + b], but may be evaluated if the segment is straight or slightly curved. 
This s-dependent value of course still is to be integrated over all s. 

To get an approximation for L in the straight wire case use cylinder coordinates with 
a length s along the axis and area element dA = rdrdcj) (see figure This leads to 



dsL{s), (A.l) 

Lo [s) = ( I , ds'dA'dA] 

\J |x(s,r,0)-x'| /|,(x')-«|<6 

In the latter integral x extends over the plane through the centre of a cylindrical segment, 
x' extends over the complete segment. The integral Lq (s) of course is independent of 
s. The integral over s' (from —b to b) may be performed using |x — x'p = A^^ + s'^, 
Ar2 = + r'2 - 2rr' cos (0 - 0'), 

Lo(0) ^ 2 j dAdA'j{r)j (r') asinh (b/N) 

= 2 j dAdA'i (r) j (r') {In {2b/ a) - In (N/a) + A^ {N/b) + ...} 

In the second line the expansion 

asinh (x) = In (2x) + (1/x) (A.2) 

A(x) = i^^4^^tlClx- = xV4 - 3x732 + ... 

IV ' L^n=\ 2-4...2n 2n ' ' 
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was used. The expansion converges because of = O (a) ^ b. It doesn't matter 
which 0' occurs in the 0-integral and thus we now set 0' = 0. Because of J dA |j| = 1 
the leading term simply becomes 21n(26/a). The second term follows from ln(A^/a) = 
\ In (p^ + p'^ — 2pp' COS0) with p = r /a and p' = r'/a and 

/ In (p2 + p'2 _ 2pp' cos 0) # = 2 In (p>) , (A.3) 

where In (p>) = 5 (p — p') Inp + (p' — p) In p'. This term thus vanishes in the skin effect 
case where the current differs from only for p = p = 1. The current density in the 
constant current case is j (r) = 1/ (vra^) and the second term becomes 

-^2 [\pp f dp'p'\n{p^) = 1/2. 



A good approximation for b ^ a thus is 

Lo(0) = (0) + O (aVfe') , (A.4) 
L^(0) = 21n(26/a) + 2r 

with Y = 1/4 for a constant current distribution and Y = in the skin effect case. The 
trivial integration over s leads to Lq = IL^ (0). 



B Semi-infinite straight segments 



One error term in equation (12) is the difference of the volume integral L{Q) and the 
curve integral (0). For current loops consisting of straight segments and sharp corners 
this difference contains the potential of straight segments of length c with a current in 
the bulk of the wire and an opposite current on the axis. These currents are collinear 
and scalar potential theory applies. This allows to generate the potential of the segment 
by subtracting the potential of a semi-infinite current distribution of this type from the 
potential of an identical semi-infinite current distribution shifted by length c along the axis. 
The essential result is that the potential of a straight segment of this type is generated by 
multipoles at the ends of the segment. 

The potential of a semi-infinite segment in cylinder coordinates is 

0(r,^) = (^[n^ + {z' - zfY^'^ - (r^ + {z' - zfY^'^"^ dz' 

= (in^ + asinh^-asinh^). 

The segment extends along 2; < 0, A^^ = + r'^ — 2rr' cos 0', and the angular bracket (...) 
is an abbreviation for j d(j)'r'dr'j (r') ... with r' < a. The current j (r') is a unit current 
(1) = 1. To verify that the potential concentrates around the end at 2; = first consider 
the case r = and z > a > r'. For r — )■ one gets N = r' and the above equation leads 



to (see equation (A.2)) 



Hr = 0,z) = {ln'-^- asinh ^) = + -) = ^ + « (aVz^) • 
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In the skin effect case of course (r'^) = a^. The potential thus decays hke l/z"^ in the 
direction of the axis. 

Next consider the cones z = qr, with q ^ some constant and r > a, 

-1 / r' \ q 

(j) (r, qr) = ( — In I IH — 5"~2 2 cos (j)' ) + asinh q — asinh 



2 V '"^ / v^l +r'2/r2 - 2 cos (00 ^'/^ 

The integral J d(f)' In (...) vanishes (this corresponds to Gauss's law), and expanding the 
asinh in 1/r leads to 

The potential thus decays in all directions like a2 |x| , and a straight segment of finite 
length corresponds to two such multipoles, one at each end. 



C Contribution from a curved segment 



The goal is to evaluate the integral L{0) from equation (11) for a segment of length 2b 
and constant curvature radius R, 

^"'"'-y |x(0.r»-x-| ''"'-|' l>^''''^''- '-°<fi 

The distance up to order O [R"'^) follows from 

(x-x')' = s'2 + Ar2 + g, 

q = s'^ [(r' cos (p' + r cos 0) i?"^ + rrR'^ cos cos 0'] - s'V (l2i?^) + s'V (360i?^) 
Expanding in g gives 

11 q 3g2 

+ r^ + - 



|X-X'| (s/2 + ]V2)l/2 2 (S'2 + Ar2)3/2 g ^^,2 + ;V2)^/2 

The in the denominators is negligible for b ^ a and the elementary integrals lead to 

£H(o)-£o(o)-ii^ + o(||.„^). (C.l) 

The term of order O {a2/ R^) contributes to the error of the self inductance formula with 
the expected order of magnitude. It is essential however, that the O {b^/R"^) term, which 
is of order a/R because of 6 = y/oR, cancels against a contribution from the curve integral 

(0) from equation (10). 

To verify this cancellation start with 

±j ^ -cos(j(0),j(6)) _ -cos(j(0),j(-&)) 
db ' |x(0)-x(6)| |x(0)-x(-6)| ■ 
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m' m o„. 




Figure 3: Two non-adjacent coplanar straight segments. The coordinates m, m' and n, 
n'measure the distance of the end points from the intersection of the segment extensions. 

Inserting \x — x'\ = 2Rsmb/2R and cos (j (0) , j (6)) = cosb/r for a section with curvature 
radius R gives 

c^^ -2 / 11 / b^^ 



with integral 



— 2b 11 b'^ 

L, (0) = const - 21n - + - — + ... (C.2) 

D Curve integral for adjacent straight segments 

For completeness we display here the curve integral contribution to the self inductance 
from adjacent straight segments of length c and d with an angle a between the currents, 

T / -, \ Mo / N r . . d + c cos a , . . c + d cos a 

L^[c,d,a) = — cos(a)|casmh haasmh — ; 

2n c sin a d sin a 

, ,\ . 1 cos a 2b . , 1 — cos a , 

— (c + a) asmh ^^^^^^^= asmh j. 

sma A/l(l-cosa) sma 

For each corner such a term is to be added to the contribution ([T]) of the (straight) segments 
by themselves. The 6-term is of order O (a) for b = a/2 and normally may be neglected. 



E Curve integral for non-adjacent coplanar straight 
segments 

What then is missing for calculating the self inductance of a loop consisting of arbitrary 
coplanar straight segments is the mutual contribution from non- adjacent straight segments, 
see figure (|3|. 

L (m, m', n, n', a) = — {A (m', n', n,a) + A (n', m', m,a) + A (m, n, n', a) + A (n, m, m', 

27r 

. , , . , f u + w COS a \ . , f V + w cos a \ 

A (w, M, f , a) = w asmh — w asmh 

\ w sin a / \ w sin a / 

This leads to an unwieldy expression already for a hexagon, but the calculation of the self 
inductance of such loops is a matter of algebra and geometry. 
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